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I. PROBLEM DESCRIPTION 


aia _ 
A : This program uses an improved linear analytical model for the calculation 
of launch vehicle steady-state response to applied sinusoidal loads acting in 
the longitudinal direction. The present approach utilizes a finite element 
technique to construct the total launch vehicle stiffness [K] and mass [M] 
characteristics by subdividing the prototype structure into a consistent set 
of (1) axisymmetric shell components a, to represent as separate units the 
fairing, interstage structure, bulkhead, tank walls and engine thrust structure; 
(2) fluid components b, and (3) mass-spring components ce, to provide the 
inertial and stiffness characteristics of the equipment and engines and vehicle 


supporting structure. 


The total vehicle characteristics are obtained by superposition of the 
stiffness and inertial characteristics of the individual shell, fluid and mass- 
spring components which are computed using a generalized coordinate approach. 
®luid force and inertial coupling between all structural components is 
accomplished by assuming fluid motions consistent with the shell component 
distortions. The superposition technique automatically assures displacement 
‘compatibility and satisfies force equilibrium at the joints between components. 
After the complete system matrix has been formulated, displacement boundary 
conditions are introduced by removing appropriate rows and columns of matrix 
coefficients corresponding to points on the vehicle and its supports which 


are rigidly restrained from motion. 


The coupled system natural frequencies and mode shapes are obtained from 


the eigenvalue equation constructed with the total stiffness and mass matrices 


Tx {o} - p> [M1 {w} = 0 


in which p is the circular frequency of the system and {a} is the natural mode 
vector whose components are the longitudinal, radial and rotational displace- 
ments at discrete points on the vehicle. The steady-state response due to 
simple harmonic loads is determined using a standard modal response procedure 
which expresses the total displacement, velocity, acceleration and stress 


responses as the linear superposition 9f the individual modal responses based 


on an assumed modal damping. ) he 4: an POZE 63203/] | 
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II. 


PROGRAMAING ASPECTS 


This section contains the following programming information: 


1. Routine description 

2. Tape format 

3. Overlay program structure 

4. Overlay core storage allocation 
5. Deck set-up 


6. Flow charts 


II-1 


Name 


MAINE 


MATC 


MATIV 


MATO 
MATL 
MAT2 
MAT3 


MAT5 


MATO 


MATT 


MATS 


MATS 


MATP 


UNOS 
UNOL 
UNO2 
BLDATA 
BLDAT1 
BLDAT2 


PRC 


ROUTINE DESCRIPTION 
Description 


Main program to generate program constants and to control 
the subprograms. 


Matrix construction - Move B matrix into A matrix starting 
at (i + 1)t% colum ane (j + 1)t2 row. 


Matrix inverse control - initialize and set up matrix for 
computing inverse. 


Matrix clear - store zeros in matrix elements. 


Matrix transpose - C = Al 


Matrix multiplication - C = AB 
Matrix addition - C = AtB 


Matrix multiplication - C = AD where D is a diagonal matrix 
stored as a colum. 


Matrix multiplication - C = apa? where Dis a diagonal 
matrix stored as a column. 


Matrix multiplication - C = ADB where D is a diagonal 
matrix stored as a column. 


Matrix multiplication - C = ADB where Dis a diagonal matrix 


stored as a column. 
Matrix inverse - C = a 


Matrix print - print a rectangular matrix or the lower 
triangular portion of a symmetric matrix. 


Tape specification - generate FILE card for tape UNITO9. 
Tape specification - generate FILE card for tape UNITO1. 
Tape specification - generate FILE card for tape UNITO2, 
Block data - generate val 

Block data - generate D1, [Wy CW], and (W] 

Block data - generate [W] 


Control program for reading input data and processing 
tables. 


II-2 


Name 
DASYS 
DASH 


DAFL 
DASM 
TABLE 
FUN2 
FINP 


SCNTL 


COMPK 
GENKM 
GENUV 


FCNTL 


GENXT 


TABLL 


DIAGM 


GENFM 


VV 


Description 


Read system input data and save polynomial matrices on Tape l. 


Read shell input data and generate shell component tables. 


Read fluid component input data and store as a column in 
cémmgn /DASH./ 


Read spring-mass input data and save stiffness matrix [K1,, 
mass matrix (M]e, and identification vector on tape 2. 


Generate shell component tables which are functions of F at 
16 points and compute initial stresses. Save processed data 
and geometry data on Tape 3 and Tape 4, respectively. 


Function subprogram to compute a part of initial stress. 


Input subroutine. See separate write-up for further 
description. 


Control program for generating stiffness matrix [Klg and mass 
matrix [M]g for shell components. These matrices are saved, 
along with the I.D. vector, on Tape 2. Transformation matrix 
{T] is saved on Tape 9. 


Compute 13 diagonal matrices which are functions of & and 
store results in C¢MM@N/WORK/ 


Generate stiffness matrix [K]lg and mass matrix [M]g for a 
shell component. 


Compute matrix [UV] 
Control program for generating mass matrix for fluid 
components. These matrices are saved, along with the I.D. 


vector, on Tape 2. 


$ an e 
Generate matrices Cel [ElLo» [5 1.39 Cela. 
or [1,03 [ai9 7 Fail (B50 = Eo)» (F 29 - 631)3 
ae Ne Meee 
Generate tables that are evaluated at 21 points of €&. 


Compute 35_ diagonal matrices that are used in generating 
matrices [M, ] and [M, ]- 


Generate matrix (™] for a fluid component. 


F =1 “a <3 =1 33 
Subroutine to compute [V1 [V1 ive (V1, or [v1 
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Wame Description 

WWW Compute the sum of three terms in matrix (M1 for case 2 or 
case 3. 

YYY Compute the sum of two terms in matrix (MJ and M1 

GENMB Compute mass matrix (M1, for a fluid component. 

TOTKM Control program to generate total stiffness matrix [K] and 


total mass matrix [M] from component matrices. [K1 and [M] 
are saved on Tape 4+ and Tape 3, respectively. 


BUILD Subroutine to build a part of a total matrix from a component 
matrix and I.D. vector. 


EGCNT Control program to compute natural frequencies and mode 
shapes. The system data required for response computation is 
saved on Tape l. 


FREQ Compute and print natural frequencies, mode shapes, velocities, 
and accelerations. Natural frequencies and mode shapes are 
saved on Tape l. 


EG2IM Subroutine to compute eigenvalues and corresponding eigen- 
vectors. See separate write-up for further description. 


RCNTL Control program to compute steady-state response. Data is 
taken from Tape 1 and Tape 2. 


RUSPL Compute and print response amplitudes and phase angles. 
STRES Compute and print component force amplitudes and phase 
angles. 
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IBMAP SUBROUTINE 
RW FINP 
Pege 1 


Identification 
RW FINP - Decimal, Octal, BCD, Variable Data Input 


Purpose 


To read a set of Hollerith punched data and/or header cards into core with 
one CALL statement. 


To convert the data fields to binary end store in core sccording to their 
associated conversion codes. 


Restrictions 
This routine uses FRWD to accomplish the BCD card image read and write. 


This routine uses FXEM error code 60 in case of errors such as non-H>llerith 
characters, data out of range, illegal format, subscripts too large for the 
array previously defined, etc. Upon detection of any error, control is 

sent to FXEM efter the end-of-case is encountered. 


Method 
Decimal numbers are converted to binary integers and then scaled to the 
indicated power of ten. 


Octal numbers are converted to binary integers. 


Hollerith words are stored directly. 

Range: Decimal to floating binary conversion 10*38 
Decimal to fixed binary; 1 to 9 digits** 35 
Decimal integer to binary integer; O to @-" -1 


Octal integer to binary integer; O to 239 -1 


Card images are optionally listed as encountered. If cards are listed, 
a page is ejected before listing the first card image. 


¥*The magnitude of the number depends upon the location of the decimal 
point. 
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1. The data card format, available on keypunch form Xl, X2, or X3, 
consists of four subfields containing the conversion code, location, 


number and exponent respectively. 


Card Columns 


Data Field Data Field Data Field Data Field 


Sub field 1 ae a a 
Conversion code “rT 19 37 55 
Location 2-6 20-24 38-he 56-60 
Value 7-16 25~34 43-52 61-70 
Exponent 17-18 35~36 53-54 71-72 

VALUE exp fia LOC. VALUE ExPh LOC. VALUE exe fui) LOC. VALUE EXP 1D 


2222212222222222|2 2121/2 2222(2222222222\2 2122222222 
3333313333333333/3 31313 333 31333333333 3/3 333333393 
4444444444444 414 Alte O44 al 44444 4d ld 484484 ga 
9955515555555 55 55 5191555551 55555555 5/5 555555555 
BEG EGESEGEGE6E G6 FBS C EGG 6G6G66666 G6 666666666 
TUTTI TTT TTT TT Ta) 
888888888888 88 8/8 818/88 88 888888 88 8B Bis 88 8 Bb 8 BBB 


Bis 888 6)8 88868888 ole ss oe8 sb 88 eee 88 als 8 
alggasaiggggsaggagsls a 


8120 21 22 23 2425 26 27 26 28 30 31 32 33 34135 36197138 39 40 41 42143 44 45 46 47 48.4950 51 52 3 sels 57 98 59 GOI 62 63 64 65 66 67 66 9 70171 72173 76 75:76 


99999)999999999 9/9 9191999991999999999 919999999999 
n 


78 79 BO 


eS 
pew?) 
ew 
~w 
ow 
=o 
aewws-) 


y 9 10 19 12 13 14 15 18147 18 
ML 22184 SPACE TECHNOLOGY LABORATORIES INC.~COMPUTATION AND DATA REDUCTION CZNTER 
\ eens) Oaes wii ak. : eee 
Data Field 1 Data Field 2 Data Field 3 Data Field 4 
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where conversion code is one of the alphabetic characters defined 
below which specifies the type of conversion to be used on the value 
field, the location specifies the cell into which the converted 
value field is to be stored, the value subfield contains the data 
to be converted, and the exponent contains the power of ten by 
which floating data is to be scaled, or the location of the 

binary point of fixed point data. 


The header card format consists of a conversion code in column 1, 
a sequence number or symbolic location in columns 2-6 and any 
Hollerith information in colums 7-7e. 


Decimal Points: 


Decimal points may be placed anywhere in the value field except that 
they may not occur in the same column as a minus sign (11 punch) since 
this results in & non-Hollerith character. If the decimal point 
would normally appear at the right of the number punched in the value 
field, then it is optional. 


Minus Signs: 


Minus signs are 11 punches over any digit of the field. If all of the 
available columns of the field are not used, minus signs may be punched 
as the left character of the field. 


Values: 


Values must always be written to the extreme left of a field. It is 
not necessary that the entire field be filled as the first blank 
denotes the end of value. Superfluous low order zeros should be 
omitted as they increase conversion error. 


‘The only exception to partial fields is BCI data where the entire 
field, including blanks, is stored. 


Location: 


The location may be specified by either absolute octal, a variable or 
array name, or the element subscripts in a one or two dimensional 
array. If the locations contain five digits, it is interpreted as 
octal. All five columns must be punched for octal. locations. 
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If the location contains at least one (1) non-numeric character, it is 
interpreted as a variable or array name which must appear exactly as 
given in the CALL statement (see Calling Sequence below). The contents 
of the number and exponent fields, if they are numeric data, are stored 
in the cell for the variable or the first cell for the array. This 
location then becomes the origin for all subscript locations following 
until another variable or array name is encountered. Caution must be 
taken to load an array name prior to subscript locations. 


If the location contains four or fewer digits, it is interpreted as a 
subscript except for conversion code H explained below. Single 
dimension array subscripts must be left justified with leading zeros 
optional. ‘Two dimension array subscripts must be denoted by two sub- 
fields containing i and j respectively. The i and j subfields must 
be separated by a comma. Each subfield may be 1 to 3 digits but the 
total field, including the comma, cannot exceed 5 characters. 


If the location is left blank, then the location counter within the 
routine is increased by 1 and the associated number is stored in the 
cell immediately following the cell where the last number was stored. 
Thus, an entire array may be read in by specifying the initial 

location only. 


Conversion Codes: 
Blenk: Floating decimal 


The number in the value field times the power of ten in the exponent 
field is converted to floating binary. Checks are made for overflow 
and format errors. 


F: Fixed decimal 


The number in the value field is converted to fixed point binary and 
stored with the binary point located at the position specified by the 
number in the exponent field. An overflow error check is made. 


i: Decimal integer 


The number in the value field is converted to a fixed point binary 
integer with the binary point following position 35. The exponent 
field is ignored. A decimal point is considered an error. 


B: Octal 
The value plus exponent fields are converted as a logical octal word. 


It is not necessary to include leading zeros but the first octal digit 
must always occupy the leftmost position of the field. 


II-8 


D: BCI Data 


The contents of the value plus exponent fields are interpreted as two 
BCI words and stored in two consecutive cells whose origin is specified 
by the location field. 


H: Heading card 


A card with an H in colwm 1 is considered a BCI heading card. If the 
location field is blank, the card is ignored. If the location field 
contains a one to four digit positive decimal integer V, (octal, or 
negative values not permitted) columns 7-72 of the card are stored 
directly in 11 consecutive words beginning at the current origin 

+11 * (V-1). If the location field is symbolic it becomes the current 
origin and is equivalent to V=l. 


The last variable or array named in the CALL statement is used to 
initialize the current origin so that symbols are not required on 
headers which are first in the deck. Each card may be used as one 
record of output using FORMAT option A with column 7 of the card 
providing the code for printer spacing. 


A: Variable names as data 


The value plus exponent fields are interpreted in a pseudo MAP instruc- 
tion format AAAAA T DDDDD P where the fields to replace are address, 
tag, decrement and prefix respectively. The address and decrement 
fields are defined normally to be 5 characters and the tag and prefix 
as one octal numeric character each. Any field containing less than 
the normal number of characters must end with a comma while fields 

of normal length must not. Any address or decrement field containing 
less than 5 numeric characters is converted as decimal while those of 
all 5 numeric characters are converted as'octal. Any address or 
decrement field containing at least one non-numeric character is 
interpreted as a variable or array name. Variable addresses cause the 
entire word from the compiler generated calling sequence to be loaded 
into the location word (i.e., the PZE X is stored in the location 
specified if X is the variable appearing in the address field). 
Variable decrements cause the right-most 18 bits from the compiler 
generated calling sequence to be loaded into the prefix and decrement 
of the location word. Numeric tags and prefixes are loaded directly 
into the corresponding parts of the location words. Null fields 
_are not loaded. The first blank indicates the end of the loading; 
the address only, the address-tag, the address-tag-decrement, or the 
entire word may be loaded as desired. 
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G: Temporary Origin 

The value in the first location field on the card is used as 2 temporary 
origin for tables. The location is saved and if data cards follow with 
blank location fields the corresponding data is stored consecutively 
peginning with the cell specified in the location in the © card. 

Columns 7-72 are ignored and may be used to identify the table. 


The first non-blank location starts a new origin. If this non-blank 
location is a subscript, it references the last variable or array 
named, which may or may not have been on the G card. 


J: Halt and Jump 


The location specified with this prefix must be an octal address and 

is the only part of the data field that is interpreted. The subroutine 
causes a transfer to the octal location specified and does not interpret 
the remaining fields on the card. 


: Two dimension array i 
L Y *nax 


The location field contains the name of the array to be loaded. The 
value field is defined to consist of 2 subfields, separated by a 
comma, containing the decimal integers for A ay and dD rivsg 


respectively where A cis and : eer generally appear in a DIMENSI@N 


statement. Until redefined the i and j values are retained 
max max 


a one definition 


to compute the successive subscripted locations. Blank address 
fields may.follow this array definition if successive elements 
of the array are to be loaded. 


: Two dimension array i j 
M Y tmax? Ima 


Conversion is identical to L except the entire array is preset to zero. 


definition 
x 


E: End case 


This defines an end-of-case and control is returned to the FORTRAN 
object program. The rest of this field and the remaining fields on the 
card are ignored. 


K: Omit loading the next N cells 


The integer N in the value field is used to increment the next storage 
location by N. 
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C: Preserve the location of the next storage cell. 


The location of the next available sequential storage cell is stored 
into the TABLE of symbolic information in the address of the word in 
the calling sequence corresponding to the word which contains the 
BCI representation of the symbol in the value field. 


R: Retain base location 


The location of the symbol in the value field is saved as a base 
location for relative addressing. The base is initialized to the 
last variable or array in the calling sequence. 


N: Preserve the relative location of the next storage cell. 


In addition to the function of the C code being performed the 
location of the next available sequential storage cell relative to 

the symbol defined by the last R code will be stored according 

to the next table entry. The next symbol in the TABLE will be checked 
for being null (ell blenks). If it is, the relative location 

(index value) will be stored in the associated location. If the 
symbol is not null an error is noted. 


Calling Sequence: 


The following two types of CALL statements may be used: 


I. CALL FINP (+ n,X,Y,ZETA,...,mHX (GS) ¥ (5) ZETA (@)...) where 


A. n is the number of variables and/or arrays in the list, excluding 
n itself. If n is negative, a card image listing is produced and 
if n is positive, no listing is provided. 


X, YyeZETA,... are the names of variables and/or arrays restricted 
to at most 5 characters each, one character of which is non-numeric. 


mis 6 times n. Hence, mH allows for 6n Hollerith characters to 
follow. 


x G) Y G) ZETA (2). -. is a list of the items previously named in 
exactly the same order with indicating the number, i, of blanks 
necessary to provide six Hollérith characters for each item. Since 
each_item name is restricted to 5 characters, the minimum value of 
of G) is (). 


CALL FINP (0) where the number of items is given as zero. This CALL 
statement must be used only after a CALL statement of type I has been 
executed. When the subroutine encounters a zero for the number of 
items, it immediately refers to the last executed CALL FINP with.a 


IT-11 


TBMAP SUBROUTINE 
RW FINP 
Page 8 


non-zero number of items for the names of the items to be loaded. 


Space Requirements 


776 cells 
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Identification 


RW EGeF Eigenvalues & Eigenvectors of Ax = \Bx 


Purpose 
To compute in single precision floating point all the eigenvalues and eigenvectors 


of the (real) system Ax = Bx, where A and B are symmetric, and B is positive 
definite. The computation of the eigenvectors is optional. 


Restrictions 


No internal checks are made for overflow or underflow. The eigenvectors are not 
normalized. The matrices A and B are destroyed by the subroutine. The size of 
the system is limited only by the amount of core storage available. 


Method 
A lower triangular matrix ve is found such that B = uM. The equation Ax = ABx 
is then equivalent to wZam™4x = AMx. We form the symmetric matrix F = wu Tam, “> 


and let y = Mx. Now we solve the standard eigenvalue problem Fy eAy using the 
Jacobi method routine NYEVV, which gives us the desired eigenvalues. The eigen- 


vectors are found by forming x = uty, Details of the decomposition ,B = uM, are 


given in Appendix A of routine RW MI2F. 


Usage 
Calling Sequence: 
Loc. Oper. Add, Tag, Decr.. 
Xx TSX EGeF, 4 
e+ 1 eee A, 0, B 
X+ 2 PZE Vv, O, N 
X<+ 3 Error Return 
o<+ 4 Normal Return 
Where: 
PZE are 
Operation of o¢+ 1 is MZE |? if eigenvectors nee wie desired. 


A is the first location of n? cells which contain the full A matrix upon 
entrance. The eigenvalues will be stored along the main diagonal upon 
exit. 


B is the first location of (n? + n)/2 cells containing the lower triangular 
portion of the B matrix. That is, C(B) = by 47C(B+1) = boy? C(B+2) = Doo 


C(B+3) = b.,,etc. 


3L 
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Usage (continued) 


V is the first location of n? cells which will contain the eigenvectors 
stored row-wise upon exit. The ith row of this matrix is the eigen- 
vector corresponding to the ith eigenvalue along the main diagonal of 
the matrix which has replaced the A matrix. 


If eigenvectors are not desired, V is the first location of n° cells 
used for temporary storage. 


N is the order of the system. 


Error return may be reached from two sources: 
a) If C(ACC) = 0, the matrix B is not positive definite. 


b) If c(ACC) ¢ 0, the matrix F is ill-conditioned under the criterion 
described in NYEVV as "restriction 1". 


Coding Information 


A square root routine (RWSQ2F) begins at EG2F + 3h5. 
A Jacobi method routine (NYEVV) begins at EG2F + 370. 


Mm? was formed by RWMI2F which begins at EG2F + 216. 


Space Requirements 


Program - 706 cells 
COMMON - (114N) cells 


Timing 
Approximately .OO5N> seconds, if eigenvectors are desired, or .0035N> 
seconds, if eigenvectors are not desired. 
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TAPE FORMAT 


Five tapes are used to store intermediate results, to transmit data 
between program links, and to save information for computing steady-state 
response. The following is a general description of the contents of each 


tape: 


Tape #1 

This tape is used as a scratch tape in subprogram DASYS to store the 
polynomial matrices [A], and (BI. which will be used in subprograms SCNTL and 
FCNTL to compute the component stiffness matrix and mass matrix. It contains 


2 x N. number of records as shown: 


P 
fA]. (BI. k= °152,3; «a5 Np 


Subprograms EGCNT and FREQ regenerate tape 1 to be saved for steady- 


state response computation. It contains four records: 


Record No. 1 - Processed response data that are stored in CéMM¢N/SYSO/: 
HHEAD, P,, C,; Nhe Nop? fis “fs, m; 


Record No. 2 - System input data that are stored in cémMgn/sys1/: 
Nw Nop Np» Ny 8, No» Ni Np» &> Ney 
Ny opt,» opt,, opt, and opt). 


Record No. 3 - Natural frequencies Py in radians per second. 


Record No. 4 - Mode shapes {o,3 
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Tape #2 


This tape contains the stiffness and mass matrices of all the individual 
components for constructing the total stiffness matrix and mass matrix and for 
computing steady-state response. It is generated by subprograms DASM, SCNTL 
and FCNTL in the following order: 


1. Spring-mass I.D. number, matrix dimension, I.D. vector, 


stiffness matrix and mass matrix (2 x Ny records): 


Cy» ns (IDC), , (KJ. (M], i= 1,2,3, Vee Nu 


2. Shell component I.D. number, matrix dimension, I.D. vector, 


stiffness matrix and mass matrix (2 x Ng records ): 


ID), (kK, 1,, (J, 


3. Fluid component I.D. number, matrix dimension, I.D. vector, 


and mass matrix (2 x Np records): 


Tape #3 


This tape is used to transmit the processed shell data from subprogram 


TABLE in link 1 to subprogram SCNTL in link 2. It contains 3 x Ng number of 


records: 
(DAS1), 
Record No. 1 - a is the shell I.D. number and ITEM, .equal to seven, is the 


number of tables in record no. 2. 
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Record No. 2 - Seven processed shell tables and three computed constants 
that are stored in C¢MM@N/DAS1/: 
g, Yr, To Y; Ng, t, and D 


k, d, and a, 


1 


Record No. 3 - Processed tables of shell orthotropic constants that are 
stored in CéMM@N/DAS2/: 
Cll, C12, C22, C33, C34 and Ch4. 


Tape 3 is utilized again in subprogram TOTKM to store the dimension 


No - No and the total mass matrix [M] as record no. 1 and record no. 2, 


respectively. The total mass matrix is saved for the computation of natural 


frequencies in subprogram EGCNT. 


Tape ff 


This tape is used to save the shell input data generated in subprogram 
TABLE to be used in subprogram SCNTL. Each record contains shell input data 


that is stored in c@itidn/DASO/. The table contains a total of Ng records. 


Tape 4 is utilized again in subprogram T@TKM to store the dimension 


N,- N, and the total stiffness matrix ([K] as record no. 1 and record no. 2, 


C 6) 
respectively. The total stiffness matrix is saved for the computation of 


natural frequencies is subprogram HEGCNT. 


Tape 79 


Tape 9 is used to store the shell component transformation matrices 
which are generated in subprogram GENKM and to be used in subprogram FCNTL. 


It contains two records per shell component as shown: 


where a is the shell component I.D. number 
U+V is the dimension of matrix [T] 


(T] is the transformation matrix. 
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OVERLAY CORE STORAGE ALLOCATION 


LINK O 


LIBRARY SUBROUTINES 


Link 3 


Origin BETA 
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DECK SET-UP 


General 
Control Cards 


Link 0 Link 1 Link 2 Link 3 Link 4 
MAINF BLDATA SCNTL FCNTL GENXI 


MATC BLDAT1L  COMPK TABL1 Data Decks 
MATIV BLDAT2  GENKM DIAGM 

MATO PROC GENUV End-of-File 
MAT1 DASYS 

MAT2 DASH 

MAT3 DAFL 

MATS DASM $ORIGIN ALPHA 


MAT6 TABLE 
MATT = FUN2 fun / 


ae $ORIGIN ALPHA 
SORIGIN ALPHA 


PORIGIN BETA 
r——_————_— 


$ORIGIN BETA 


SORIGIN ALPHA 
$ORIGIN ALPHA i 


Link 5 Link 6 Link 7 Link 8 
$ORIGIN ALPHA 


GENIM TOTKM EGCNT RCNTL 


r VV BUILD FREQ RESP1 
pried www *EGQIM — STRES 
YYY . 
GENMB 
$IBJOB 
$EXECUTE IBJOB 

Tape Cards 

General *Source language of the subroutine 

Control Cards is MAP 
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MAINF MAIN CONTROL PROGRAM 


BEGIN 


Read input cards (including stiffness 
matrix [K]c and mass matrix [Mle of 

spring-mass components), and process 
geometry data CALL PROC 


Construct stiffness matrix 
[Ka and mass matrix [M]a 
for N, number of shell 
components CALL SCNTL 


Compute and print steady-state 
response, velocities, accelera- 
tions, and component forces 


CALL RCNTL 
alta 


Construct mass matrix [M7b 
for Np number of fluid 
components when Np # 0 


CALL FCNTL 


Is this 
the last 
case? 


Construct total stiffness 
matrix [K] and total mass 


matrix [M] from component | | sT¢P 
matrices qary, TOTKM 


Compute and print 
frequencies, mode shapes, 
velocities, and accelera- 


tions ary, EGCNT 
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PROC CONTROL PROGRAM FOR DATA INPUT AND PROCESSING 


a 


Read system input data and save 
polynomial matrices [A}, [Blk 
on tape 1 

CALL Dh ee cee ee eee 


Read shell input ae input data and save 
on tape 4. Generate shell 
component tables and save on 
tape 3. CALL_DASH 


Read fluid component input data 
and store in 
COMM@N DAS4 

CALL DAFL 


Read spring-mass input data and 

save stiffness matrix [K ]e and 

mass matrix [M]e on tape 2. 
CALL DASM 
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READ SYSTEM INPUT DATA 
DASYS AND SAVE POLYNOMIAL MATRICES ON TAPE 1 


START 


Clear input regions 
and rewind tape 1 


Read system input data 
CALL FINP 


Generate a complete vector 
for modal damping factor 
if Nyp 4 0 


S>=0 
RETURN 


Read a pair of polynomial 
matrices [A] and [B] 


CALL FINP 


Save matrices [A] and 
{B] on tape 1 


Last 
pair of [A] 
and [B]? 


no 


Rewind tape 1 
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DASH 


READ SHELL INPUT DATA AND 
GENERATE SHELL COMPONENT TABLES 


Rewind Tape 3 and 4 


Read shell input data 
CALL FINP 


Generate shell component tables 
and save on tape 3. Save shell 
input data on tape 4. 


CALL TABLE 


Store shell component 
I.D. vector in 


COMMGN DAS3 


How has this been done 
for all shell components? 
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DAFL 


READ FLUID COMPONENT INPUT DATA 


(| START } 


Is there 
any fluid 
component? 


Read fluid component 
input data 


CALL FINP 


Store input data as a 
column in 


COMMgN DASL 


' Has this been done for 
all fluid components? 
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DASM READ SPRING-MASS INPUT DATA 


START 


Is there 
ey spring-mass 
component? 


Read spring-mass input 
data 


CALL FINP 


Save stiffness matrix [K],, 
mass matrix [M],, and I.D. 
vector on tape 2 


Has this been done for 
all spring-mass components? 


RETURN <j 
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TABLE GENERATE SHELL COMPONENT TABLES 


START 


Generate tables for othotropic 
shell constants C1, Cia, Coo, t, 


C33, C325 and Chh 


conical shell a 20 a <0 ellipsoidal shell 


| Set k = 0 | Compute a,b, and k 


Set Cl=o-" 
C2=Go 


Generate tables 
@= Cle + C2, 
Fy »Yp»T, and D = 1/(r, Ch) 


Generate tables 9, rj, 
ro, r and D 


Compute d, and qd, 


1 


Compute qa, and a, 


Compute initial stress 


3 Compute initial stress 
N 


Nd 


Write generated tables (C@MM@N DASL) 
on Tape 3 


trite input data (C¢MM¢N DASO) 
on Tape 4 
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SCNTL CONSTRUCT STIFFNESS MATRIX AND MASS MATRIX 
FOR SHELL COMPONENTS 


Rewind Tape 3, 4 & 9 
Rewind Tape 1 


For the ith shell, read in tables from 
tape 3 and input data from tape 4. 


Search and read the corresponding poly- 
nomial matrices [A] and [B] from Tape 1 


Compute 13 diagonal matrices 
[KM]; 6 x 16 M=1,2,----,13 CALL COMPK 


Generate shell component stiffness 
matrix and mass matrix CALL GENKM 


Write shell component stiffness matrix 
and mass matrix on Tape 2. 


Has this been 

done for all 
shell 

components? 
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GENKM GENERATE STIFFNESS MATRIX AND MASS MATRIX 
FOR A SHELL COMPONENT 


Compute U1, [V,1, fo Ww (U,] ana FV,1 


! 
| Compute ik, re | 


i 


Y 
Compute [K,,1 = (K.,] * Kept 


Compute iM, Pe 


Compute fUV? 
CALL GENUV 


Compute local transformation matrix [T] 


RETURN 


Compute Key] = rye [k,,] rr 


Compute (M1 = rr qe m1 (T] 


. a . 
Compute scaling factor — for mass matrix 


Me, 


Print shell. stiffness 
matrix with equili- 
brium check and shell 
mass matrix with mass 
check 


Is 
Input Mass 
Zero? 


Apply ceatine factor to mass matrix 
Mia = 2 ne 
Ma 
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FCNTL CONSTRUCT MASS MATRIX [M}. FOR FLUID COMPONENTS 


START 


Pick up input data of the 


ith fluid elements from 
storage 


Search and read in input data of 
associated shell components from 
Tape 4. 


Search and read in associated poly- 
nomial matrices [Aly (B)..> [Alo 


[Blo [A].3 and (Bla from Tape l. 
Should fluid component 
mass matrix be zero? 
Set MI, = 0 


Generate matrices [€]al, [e]a2, 
[e]e3, [£]e2 or [£]a23 (E19 - 844), 


~ e op . 
(E59 Bods (B30 53)3 Vy Yo? 


aN CALL GENXI 


Generate tables that are evaluated 


at 21 points of F (m=1,2, 3} 
CALL TABLL 
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©) 


Compute 35 diagonal matrices that 
are used in generating matrices 
(My J] and [Mo] 

CALL DIAGM 


ompute matrix [M] 
CALL GENMF 


Search and read in associated 
transformation matrices T1, To, & 


» from Tape 9 


Compute fluid component matrix 
(M4, CALL GENMB 


Compute mass check if input fluid 
mass is not zero and condition mass 
matrix for symmet 


| Print fluid component matrix 
[M}p when I@PT1 = 1 


Write fluid component matrix 
[M}, on Tape 2. 
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GENFM 


COMPUTE MATRIX [M] FOR A FLUID COMPONENT 


ny 


Compute [Up], [Vey ]> Mol (Vypls Czy] and [¥,5] 


ol) [V0] for case 3 


| Compute (Uo 1, (Vv ] for case 2, or (U,. 


Clear [M] and compute (0,1) 


Compute (v1) 


teen 


| compute (¥,17 if this is case 2 


Compute the lst part of (mM, + [M1 and add to (M1 


Compute [u,,2 7 


Compute [V2] 


+ (M1 and add to [M] 


Compute the 2nd part of [M, ] 
al 


| Compute (0,31 | 
| Compute (V,31 


Compute C¥,3] if this is case 3 


Compute the 3rd part of (Mi, 7 + (M1 and add to [M] : 
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GENMB 


COMPUTE MASS MATRIX FOR A FLUID COMPONENT 


START 


Assemble local coordinate transformation 
matrix [T] from the associated shell 
component transformation matrices [T,], 
[To1 and [T31 


Compute [M]} = cor (M} (PJ 


Is 
input mass 
zero? 


yes 


no 


Compute scaling factor 


mass matrix. 


Apply scaling factor to mass matrix 


M 
tM} = 2 om 
hip 
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TOTKM CONSTRUCT TOTAL STIFFNESS MATRIX 
AND TOTAL MASS MATRIX FROM COMPONENT MATRICES 


Rewind Tape 2,3, & 4 


Clear total stiffness 


matrix & total mass 
matrix. 


Read component stiffness matrix, 
mass matrix & I.D. vector from 
Tape 2 


Expand the component stiffness 


matrix & mass matrix, & superim- 
pose on the total stiffness & mass 


matrices. CALL BUILD 


Has this been done for 


all spring mass & shell 
components? 


Read fluid component mass matrix 
and I.D. vector from Tape 2. 


“xpand the fluid component mass 
matrix and superimpose on the tota 
mass matrix CALL BUILD 


© 


Write total stiffness matrix and 
total mass matrix on Tape 4 and 
Tape 3, respectively 


Print the total stiffness matrix & 
the total mass matrix when 
T@PT2 = 1 


4 Rewind Tape 2,3,& 4 
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CONTROL PROGRAM TO 
EGCNT COMPUTE AND PRINT FREQUENCIES, MODE | SHAPES, VELOCITIES, 


AND ACCELERATIONS 


Rewind Tapes 1,3, & 4 


Read total stiffness matrix 
[K] from Tape 4 


j 


Read total mass matrix [mM] 
from Tape 3 


Compute eigenvalues and 
eigenvectors CALL EG2kM 


Save input data on Tape 1 for 


steady-state response 
| computations. 


| Compute and print frequencies 
| bo-3 -p-s.), mode shapes, 
velocities, and accelerations 


| CALL FREQ | 


le Sak 


Rewind Tapes 3 & 4 | 
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FREQ 


COMPUTE AND PRINT FREQUENCIES, MODE SHAPES, VELOCITIES, 


AND ACCELERATIONS 


START 


Rearrange eigenvalues We and 
corresponding eigenvectors {ii} 


in ascending order 


Compute and print frequencies 
in cycles per second Pe wy, /2n 


Compute and print mode shapes 


{osx} 


Compute and print mode 
velocities {“ik} wy, 


Compute and print mode 


accelerations -{%ik} we 


Save frequencies We in radians 


per second on Tape 1 


RCNIL COMPUTE AND PRINT STEADY-STATE RESPONSE, VELOCITIES 
ACCELERATIONS, AND COMPONENT FORCES 


Rewind Tape L 


Read system data and 
response data from tape 1 


Read frequencies (Rad./Sec.} and 
mode shapes from tape 1 


Generate load vector from input 
loads and I.D. vector 


Compute the response amplitudes 
and phase angles 
CALL RESPL 


Compute modal velocities and 
modal accelerations 


Print response amplitudes and phase angles, 
modal velocities, and modal accelerations 


Is computation of 
(res) component forces required? 


Compute and print the component 
forces 


CALL STRES 
Has this been done for 
(aio) all input frequencies? 


(es) Bs RETURN 


Compute the product of the mode 
shape and the applied load 


Q, = {o4,) {P} 


Compute 
2 2 


2,2, 2 2, 2.5 


Compute Py, 
5 = tan 27k (zw) | 


Ses 


Has this been done for 
all the modes (k =n )? 


Compute the steady-state response 
ainplitude Ry 


Compute the steady-state response 
phase angle 04 


Has this been done for 
all the points (i = n)? 


( RETURN } 
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STRES COMPUTE AND PRINT COMPONENT FORCES 


Rewind Tape 2 


Skip spring-mass stiffness matrix 
(K], and mass matrix [M], on 
Tape 2 


Read shell component stiffness 


matrix (K], from Tape 2 es aa 


Compute force components 


x, = Kos), {R, sin 55} 


Y; = KJ, {R, cos 63 


Compute force amplitude 


2 7 as 
poe 


Compute force phase angle 
5. = tan 
J 


“1X; 
Yj 


Print force amplitude and phase 
angle 


Has this been done for RETURN 
(10) all shell components ? S c ( _ReruRN 
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Tift. OPERATING INSTRUCTIONS 


A. HARDWARE REQUIREMENTS 


1. IBM 7090 or 7094 computer with 32K core. 
2. Data channels and tape units for the software. 
3. No on-line printer and punch. 


4, Peripheral equipment for card-to-tape, tape-to-printer, and 
tape-to-punch as required by standard IBM system. 


5. It is assumed that standard keypunch and verifier machines, 
card readers, and printers are used; therefore, these have 
not been specified. 


B. SOFTWARE REQUIREMENTS 
1. IBSYS Operating System Tape Version 9 


Basic Monitor (IBSYS} Version 4 

. Processor (IBJ@B) Version 3 
Assembly Language (MAP) Version 3 
. Assembly Language (FTC) Version 3 
. Loader (IBLDR) Version 3 


oa0o0n 


2. The Fortran IV 1/0 Library Subroutines 
C. TAPE USAGE 


The standard Fortran input/output configuration, as described in 
IBM 7090/7094 IBSYS Operating System: IBJOB Processor, Form C28-6275, is 
used for the logical tape unit designation and the file specifications are 


as follows: 


UNITOL FILE ,A(1), HOLD, INOUT, BLK = 256, BIN 

UNITO2 FILE ,B(1), HOLD, INOUT, BLK = 256, BIN 

UNITO3 FILE ,UT3, READY, INOUT, BLK = 256, BIN, NOLIST 

UNITOL FILE ,UT4, READY, INOUT, BLK = 256, BIN, NOLIST 

UNITO5 FILE ,IN, READY, INPUT, BLK = 14, MULTIREEL, BCD, NOLIST 
UNITO6 FILE ,OU, READY, OUTPUT, BLK = 110, MULTIREEL, BCD, NOLIST 
UNITO7 FILE ,PP, READY, OUTPUT, BLK = 28, MULTIREEL, BIN, NOLIST 
UNITO9 FILE ,, READY, INOUT, BIK = 256, BIN, NOLIST 


The logical tape unit 5 and 6 are used as input and output, 
respectively. The logical tape unit 1, 2, 3, 4, and 9 are used by this 


program as intermediate storage tapes and tape SYSLB2 is assigned to overlay. 
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D. DECK SET-UP 


1. Installation control cards (e.g. $DATE card). 


o. $EXECUTE IBJGB 
$IBIGB (Options) 


Program 
decks 


: 7,8 punches in Col. 1 & ) 
End of File 12,1,4,7 " ” tt 3 & h 


Data 
decks 


| 


End of File ( Same as above ) 


For a description of these control cards see: 


(IBM 7090/7094 Operation Systems Basic Monitor (IBSYS), 
Form C28-6248) 


(IBM 7090/7094 IBSYS Operating System IBJ@B Processor, 
Form C28-6275) 


3. Installation control cards (e.g. $ENDFILE) 


E. TIMING 
The computer time for each case depends upon the number of individual 
components, the total number of system coordinates selected for the vehicle, and 
the number of forcing frequencies for the response. A typical case of twenty-one 
shells, four fluids, eight spring-mass components, and seventy-eight system 
coordinates took approximately six minutes to compute the natural frequencies 
(S <0). For the same case, it took approximately six minutes to compute the 


response (S = 2, opt), = 1) for fifty forcing frequencies. 


F. MULTIPLE CASES 


The cases may be stacked for one computer run if any one of the 


following conditions are met: 
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1. All the cases have the same steady-state response option. 


2. 


ach case has the steady-state response option of either 
S8<0OorS:0. he necessary data for the steady-state 
computation is saved on tape 1 and tape 2 for the last 


case only. 


3. Each case has the steady-state response option of either 
S =lorS§=2. All cases must be using the same set of 


data on tape 1 and tape 2 for the steady-state computation. 


TT3 


IV. INPUT DESCRIPTION 


The input data is read in with symbolic data input subroutine FINP. See 


the FINP subroutine write-up for further details. 


The following is a general description of the input data. The data 
consists of four logical blocks and must be set up in the same order that is 
given on the program input sheet. The symbols used in this description are 


the same as those stated in the left-hand column of the program input sheet. 
A. THE SYSTEM INPUT DATA 
1. Heading 


HHEAD ' is one line of BCD characters which will be 
printed as the title of the printed output. 
The number of BCD words (6 characters per 


word) must not exceed 11. 
2. Input Parameters 


No is the total number of system coordinates which 


include the fixed coordinates No: 


Ny - No must not exceed 80. 
N is the total number of shell components. Ng 


cannot be zero and must not exceed ho. 


is the total number of fluid components. N 


F F 
must not exceed 6. 

Nu is the total number of spring-mass components. 
Nu must not exceed 30. 

No is the total number of fixed coordinates. 
No - Ny must not exceed 80. 

3. Applied Loads and Forcing Frequencies 

Ny is the total number of discrete applied loads. 
N,, must not exceed 80. 

C; is the applied load coordinates 


as OS Can et p Cea 
Le ree Na Ny. 
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PS is the discrete applied loads 
Pye Py» Py witereg Py . 
L 
Ny is the number of sets of forcing function 
frequencies. 
ts Af, m, f is the frequency of the forcing function 


in eycles per second. This program will 
compute the steady-state response for the 
frequencies Ps. f, + Af,» f. + OAPs are, 
.. 2 (m, - 1) Af,, 1 =1, 2, ..., M. 


Modal Damping Factors 


Nyy is the number of input n. Program will 
generate a complete table of n by setting 
n a\y prety WH equal to 

(Nant) (Nyent2) (No No) 

"Wen 

Mh is the ratio of the assumed damping to the 
critical damping in mode k 
Ths No» Ng» seas My . 

7 Daly 

Ratio of Accelerations 

gz is the ratio of the vehicle acceleration to 
the acceleration of gravity. 

Steady-State Response Option 

s is a fixed point word which controls the 


option of computing steady-state response. 

S<0 indicates that the computation of the steady- 
state response is not included. The necessary 
data for the steady-state computation is saved 


on Tape 1 and 2. 


S=0 indicates that the computation of the steady- 


state response is included. 


S>0 


7. Print Options 


opt, 


opt, 


opt 


is the option to compute the steady-state 
response only. The necessary data should be 


available on Tape 1 and 2. 


S=1, Heading (Item A.1.), Input Parameters 
(Item A.2.), and Modal Damping Factors 
(Item A.4.) are supplied as input. 


S=2, Heading (Item A.1.), Input Parameters 
(Item A.2.), Applied Loads and Forcing 
Frequencies (Item A.3.), Modal Damping 
Factors (Item A.4.), and the option 


word opt), are supplied as input. 


is a option word which controls the output 
of stiffness matrix and mass matrix of the 


shell and the fluid components. 


1, print the component matrices 


oO 
3 
ct 

it 


= 0, suppress the printing of component 


fe) 

Le} 

ct 
I 


matrices. 


is a option word which controls the printing 


of total stiffness matrix and total mass 


matrix. 

opt, = 1, print the total stiffness and mass 
matrices. 

opt, = 0, suppress tne printing of total 


stiffness and mass matrices. 


is a option word which sets the rigid body 


frequency to zero for computing the response. 


= 1, set the first frequency to zero. 


° 

Le] 
ct 
1 


) 
opt. = 0, do not set the first frequency to 
Zero. 
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opt), 


EI 


B. THE SHELL COMPONENT 
1. %I.D. Number 


a 


is a option word which controls the computation 
and printing of forces for the steady-state 


response. 
opt), = 1, compute and print the forces. 
opt), = 0, do not compute and print the forces. 


is the number of frequencies, mode shapes, 
velocities and accelerations that will be 


printed as the final output. 


8. Polynomial Matrices 


is the total number of polynomial matrices. 


is the number of rows of polynomial matrix 
[ANG x 11° 


is the number of rows of polynomial matrix 


is UL x 11 polynomial matrix. 


is ¥,. x 11 polynomial matrix. 


k= 1, 2, 3, veeey Np 


The input sequence of the polynomial matrices 
establishes the identification number k which 
is referred by the shell components. The 

subscript k is used as the polynomial matrix 


identification number by the shell components. 


INPUT DATA 


is the identification number for shell 


component a where 0 <a < Ng 


+ a indicates a conical shell component 


- a indicates an ellipsoidal shell component 
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2. 


Coordinates 


Coordinate I.D. 


(1D), 


are the total number of system coordinates. 


are the total number of local coordinates 


U,; V must not exceed 11. 
Vector 


is the identification vector which is used to 
position the elements for building total 
stiffness and mass matrices. The length of 
the vector must be equal to U + V and the 


number must not be greater than Nas 


4. Polynomial Matrix Identification Number 


k 


is the polynomial matrix identification number 
which refers to polynomial matrices [A], and 
(B},. in the system input data. 


Data 


Shell Geometric 


# 


is the meridional angle for conical shell and 
is the edge meridional angle for ellipsoidal 


shell. $, is input in degrees. 


is the height of conical shell 
+ L indicates converging upward 
- L indicates converging downward 


L = 0 for ellipsoidal shell input 


is the lower radius of conical shell 


Ry = O for ellipsoidal shell input 


is the height of ellipsoidal shell 
+ b indicates convex upward 
- b indicates convex downward 


dD 


il 


O for conical shell input 


is the radius of the base of ellipsoidal 


shell. a = 0 for conical shell input. 


6. 


(C11), 


(C33), 
(Co) 
(Cy), 


(+), 


Orthotropic Shell Constants and Thickness 


are orthotropic shell constants at two points 


p=1,2. 


are orthotropic shell constants at four points 


p = 1,2,3,4. 


are shell thickness at two points 
p=1,2. 


Mass Density and Total Mass 


is the mass density of the shell component. 
is the total mass of the shell component. 


M, # O, the ratio of the total mass M, to the 
computed mass M, will be used as the scaling 
factor for the mass matrix. When opt, = 1, 
the scaling factor will be printed as the 


mass check of the mass matrix. 


M, = 0, no scaling factor will be used for 
the mass matrix. When opt, = 1, the computed 


mass will be printed. 


Initial Stress Data 


is the depth of interior fluid. 

is the weight density of interior fluid 
is the uniform interior pressure. 

is the depth of exterior fluid. 


is the weight density of exterior fluid. 


Pa is the uniform exterior pressure 
W is the reactive force at upper edge of conical 
shell. 


+ W produces tensile stresses. 
- W produces compressive stresses 


W = 0 for ellipsoidal shells. 


C. THE FLUID COMPONENT INPUT DATA 
1. %I.D. Number 


b is the identification number for fluid component 


b where 0 <b < Np. 


2. Associated Shell Components 


As Ap, are the identification numbers of the 


associated shell components 
and a 


3 
3. Fluid Data 
H is the depth of fluid component. 
is the mass density of fluid component. 


M is the mass of fluid component. 
M #0, the ratio of the total mass M to the 
computed mass M will be used as the scaling 
factor for the mass matrix. When opt, = 1, 
the scaling factor will be printed as the 


mass check of the mass matrix. 


M = 0, no scaling factor will be used for 
the mass matrix. When opty = 1, the computed 


mass will be printed. 
D. THE SPRING-MASS COMPONENT INPUT DATA 
1. %I.D. Number 
c is the identification number of spring-mass 


component c where O <c < Ne 


n 


3. Coordinate I.D. 


i, 


consistent. 


e. Stiffness and Mass Matrices 


is the order of the spring-mass component 


n must not exceed 10. 
is n xn stiffness matrix of spring-mass component. 


is n xn mass matrix of spring-mass component. 


Vector 


is the identification vector which is used to 
position the elements for building the total 
stiffness and mass matrices. The length of the 
vector must be equal to n and the number must 


not be greater than No: 


Note: Care must be taken so that the units of the input data are 


ex 
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THE SHELL COMPONENT INPUT DATA 


SPACE TECHNOLOGY LABORATORIES, INC. 
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THE FLUID COMPONENT INPUT DATA 
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THE SPRING-MASS COMPONENT INPUT DATA 
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STL FORM toot 


V. OUTPUT DESCRIPTION 


The output of this program consists of four parts. The following is a 


general description: 
A. INDIVIDUAL MATRICES 


When opt, = 1, the following component matrices will be printed: 


1 
1. Stiffness matrix of shell component with equilibrium check. 
2. Mass matrix of shell component with mass check. 
3. Mass matrix of fluid component with mass check. 

B. TOTAL MATRICES 


When opt, = 1, the following total matrices will be printed: 


2 
1. Total stiffness matrix. 


2. Total mass matrix. 
C. MODAL SUMMARY DATA 


The following modal summary data are printed. The number of modes 


printed depends on the input parameter Nate 

1. Natural frequencies in cycles per second. 

2. Mode shapes by columns. 

3. Mode velocities by columns. 

4. Mode accelerations by columns. 

D. RESPONSE SUMMARY DATA 

The following data are printed for the modal response expression 
{R, } 7 {R, } sin @wt - {65 }) 

1. Amplitudes {R, } 

2. Phase angles {8 } 

3. Modal velocities {R, } w 


4, Modal accelerations - {R, wo 
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When opt), = 1, the following data are printed for the component force 


expression: {s,} = {8,} sin (wt - (5,3) 


5. Amplitudes (8, } 


6. Phase angles 5, } 


VI. TEST CASE 


A one-stage launch vehicle is used as a test case to illustrate the input 


data, requirements and to present the results in the final output format. 


In this example the vehicle is represented by eleven shell components, 


two fluid components and four spring-mass components. 


The input sheets and the computer output sheets are included in this 


section. 
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